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The infrared dynamics of a minimally coupled scalar field in de Sitter spacetime can be described as 
Brownian motion of a particle in a medium of de Sitter temperature Tds = The system obeys 
a fluctuation-dissipation relation and its equilibrium distribution is Maxwell-Boltzmann, implying 
kinetic and potential energies of comparable magnitudes. The transition to equilibrium is a semi- 
classical process beyond the scope of perturbation theory for interacting fields. The stochastic kinetic 
energy of the field causes de Sitter spacetime to cool down slowly with a corresponding decrease of 
the effective vacuum energy. 
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INTRODUCTION 

It has been known for a long time that spacetimes with horizons can be assigned a temperature and exhibit 
properties similar to thermodynamical systems [l[ • This connection has been extensively explored in the case of black 
holes which admit a description in terms of corresponding thermodynamical laws. Horizons also appear in cosmology. 
■ A prominent example is de Sitter space with horizon radius R = 1 / H and temperature 

One physical manifestation of this temperature is the thermal spectrum of excitations exhibited by an Unruh detector 
\ coupled to a scalar field in de Sitter, a phenomenon experienced by an observer confined within the de Sitter horizon. 
On the other hand, fluctuations of the scalar field itself on superhorizon scales, highly relevant for inflation Q, do 
l— ~~ 'I not exhibit obvious thermal properties. In fact, inflationary fluctuations have a scale invariant and not a thermal 
I ■ spectrum so there is no apparent link to thermodynamics and the temperature ([1} in this case. 
| In this letter we establish a connection between scalar field fluctuations on superhorizon scales in de Sitter and 
0^ ' thermodynamical considerations. We will be using the cosmologically relevant fiat slicing coordinates in which the 
CNI metric reads 

(N 

O ■ ds 2 = -dt + e 2Ht dx 2 , (2) 

iri . 

and assume that the the de Sitter phase starts at some definite time labeled by t = which will be assumed sufficiently 
' earlier than all times considered here. The dynamics of a test scalar field is governed by the action 
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>-i^--v(.), (3) 



^ , leading to the classical Klein-Gordon equation 

a ■ 

BV 

(f) + 3H<f ) -e- 2Ht V 2 <f ) + — = 0, (4) 

0(f> 

where a = e Ht is the scale factor. We will consider a light scalar for which the mass m <§i H. 

On long wavelengths k > aH the gradient term in (|4]) becomes irrelevant and the motion of the scalar field is over- 
damped with the the expansion rate H providing a friction term. Assuming an analogy can be drawn, let us recall 
some facts about known dissipative systems which exhibit friction on macroscopic scales. In such systems dissipation 
is generically accompanied by fluctuations with a firm relation between the fluctuation and the dissipation parameters. 
The earliest and most famous example is the phenomenon of Brownian motion described by the Langevin equation 

dv —, . 

m—+fv = !;(t) (5) 

where v is the velocity of a Brownian particle with mass m and 7 the friction coefficient. is a stochastic force due 
to the thermal motion of the medium in which the particle is immersed and whose correlator can be inferred from 
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thermodynamic equilibrium 

Mt)Zj(t'))=^T6(t-t')6 ij , (6) 

where T is the temperature of the medium and Boltzmann's constant is set to unity (fee = !)■ Another famous 
example involves Johnson-Nyquist noise in electronic circuits. If R is the resistance, L the inductance and C the 
capacitance of the circuit, the dynamics of charge Q in the circuit is governed by the equation 

d?Q dQ Q 

where t/ cx t is an externally applied voltage and is an intrinsic omnipresent fluctuating voltage related to ther- 

modynamical fluctuations of the electromagnetic fields. Statistical physics again predicts that 

(Zumu(t'))=2RTS(t-t') (8) 

In both cases the magnitude of the friction term in the dynamical equations (dissipation) determines along with the 
temperature the amplitude of the thermally induced random force (fluctuation). This is a facet of the (classical) 
fluctuation-dissipation theorem of statistical physics. 

Given the omnipresence of fluctuation-dissipation relations in nature, one can conjecture by analogy that a fluctua- 
tion term should necessarily accompany the scalar field dynamics (Q} on long wavelengths where the friction dominates. 
We thus expect that the true dynamics of the scalar field on physical scales r > 1/H (r = e Ht x) to be governed by 

4> + 3H<i>+^=m (9) 
where the amplitude of the fluctuation term should conform to a fluctuation dissipation relation 

3 \HJ 

The numerator directly mirrors the known fluctuation-dissipation relations described above. The inverse volume 
factor is required for dimensional reasons; we will find below that it is simply the physical volume of the de Sitter 
horizon i?DS = 

The above equation shows that the "soft" part of the field <fi, comprised of modes with physical wavelengths 
k/a < I/i?DS ; executes Brownian motion in an environment of temperature Tds- It should be stressed that no 
external thermal environment is assumed and (f> is in its ground state, taken to be the Bunch-Davis vacuum. The 
source of the fluctuations are the vacuum fluctuations of 4> itself, which are amplified on long scales, and the thermal 
behavior of the soft modes arises solely from the effects of spacetime curvature on the quantum field. In what follows 
we verify © and (fTU)) as the correct long wavelength description by explicitly integrating out the "hard" modes of 
the field with physical wavelengths shorter than i?os- 

The idea that the dynamics of long wavelength light fields in (quasi-)de Sitter can be described by a stochastic 
equation was first explicitly proposed by Starobinsky [3| (see also [4|]). The treatment presented here, based on out-of- 
equilibrium methods of statistical physics, clarifies and extends Starobinsky's original stochastic inflation approach. 
The agreement of the stochastic approach with Quantum Field Theory at the perturbation level was first pointed 
out in Here, we exhibit this link in full generality and at the non-perturbative level for a scalar field. We show 
that the thermodynamic analogy alluded to above is exact and the equilibrium distribution of the scalar field is 
Maxwell-Boltzmann. The latter implies an average kinetic energy comparable to the average potential energy with 
cquipartition for free fields. Such a relation is characteristic of thermodynamic equilibrium no matter how strong the 
dissipation. These results allow for an exact thermodynamic interpretation of the de Sitter temperature for the case 
of cosmologically relevant scalar field fluctuations. 



STOCHASTIC DYNAMICS IN DESITTER 



To prove the assertions made above we proceed to obtain an effective theory for the long wavelength modes of the 
scalar field. We split cj> into a long wavelength part </>(k < eaH) and a short wavelength part ip(k > eaH): $ = 4> + <p, 
where e is a small parameter that controls the separation between the long and short sectors. As long as e <C 1 our 
final result will be independent of e. 
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The split is achieved by a window function Wu, (t) in k-space to be discussed shortly. The contributions to the action 
from (j) and tp are 



S = S[<f>]+ / d 4 xa 



(11) 



where O = 



ZHdt - K 



S[4>] is the complete action for the long wavelength fields, including all non- 
linear terms involving the long field <fr. We have collected all other interaction terms in S m t- The second term is the 
free action for the short field and a bilinear long-short coupling controlled by Wkif) + 3HWk(t); bilinear coupling 
terms that do not involve time derivatives of the window function vanish. This coupling arises solely from the time 
dependence of the long-short split and describes the influence of short wavelength modes that cross into the long 
wavelength sector due to the cosmic expansion. The form of this coupling term suggests that we define the long 
wavelength sector by the window function 



W k (t) = 1 - 
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for which Wk(t) + 3HWk(t) = S(t — jj In (fc/e£f)). Thus, the coupling is effective for any mode k when k = eaH. 

The effective long wavelength action A[cj>] is defined by integrating out the short field ip in the Closed-Time-Path 
(CTP) formalism 
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«/ c 4*ftf] Jfc a 3 [yd v +^{6<t>)} 
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where the subscript C indicates integration along the closed Keldysh (forward-backward) time contour - see @ for 
an early treatment at leading order. The tp integration can be performed with the term containing £j nt handled 
perturbatively. Its treatment will be presented elsewhere. Ignoring £ mt for now we obtain 
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where x = xeiJ and we defined 



d 4 xd 4 x' a 3 (i)a 3 (^') ^(x)AT(x, x)ip{x) 



(14) 

For e < 1, 



, vr+W) v = Mix,?) = s i°fcf^ (t 

ie considering wavelengths much larger than 1/H, the relevant terms are those linear and quadratic in the field 
-0 without the gradient term. Approximating Af(x,x') — > S(t — t') d( x ~ x ) W e therefore find that long wavelength 
correlation functions can be computed using 1 



(Q(cfi)) = II/ WW] QW e lfdta3 i^ +3H ^ + ^ +i ^ . 



(15) 



Different spatial points are uncorrelated and in what follows we will be focusing on the functional for a single point. 
It should of course be remembered that in reality each "point" corresponds to a region of physical size R ~ 1/H. 
The functional integral (|15|) can be used for pcrturbative calculations. Writing 
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we obtain the free correlation functions 

( m 



/ WW)) WW)) \ _ ( Fits') 

"))) ~ \ ^G A (t,t>) 
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a?\il}<jf 



(16) 



(17) 



1 This is known as an MSRJD functional after Martin, Siggia, Rose, Janssen and deDominicis 0, 01 
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where G^ A \t,t') are the retarded and advanced Green functions 

G R (t,t')=G A (t',t) = ^^7) (e-w(*-*0 _ e -3H( t -tO^) Q(t _ t/) (18) 

and 



Thus 



+00 

W) = -^2 J dra ^ T ) G^t^^t') (19) 




(0(^(0) --^e-S^I*-*'!. (20) 



In the massless limit G R {t,t') -> 3H ^, tl) (l - e _3H (* _t ')) 6(t - i') and the variance grows linearly with time 



3Ha 3 (t') 

H 3 



\t)) m=0 ~ ^t. (21) 



It is easy to obtain the higher order loop corrections at the coincident limit. From the propagators in (I17[) and the 
vertex in (IT5|) we have at one loop 

mf) = ^(i, t) + i(-t)A / dr a 3 (r) (>(<, r) + F(r, t)) G fl (i, r)F(r, r) , (22) 

which implies that in the massless case perturbation theory breaks down when Ht ~ -^=, while in the massive case 

2 

this occurs at Ht ~ jjj?- From these results it is clear that perturbation theory is inadequate for describing the 
evolution of the interacting system for arbitrarily large times and in particular the eventual transition to equilibrium 
which we describe below. 

THE EQUILIBRIUM DISTRIBUTION 

The long wavelength dynamics represented by (IT51) corresponds to a stochastic evolution of the field <f>. This is 
made clear by performing a Hubbard-Stratonovich transformation 

e-/-^=/ {DC}e- f dt a3 (23) 

leading to 

(2W0>= J [DC\e' S dta ^ e ] J[D<p]Q(<t>)5U + 3H<i>+^-^ (24) 



Thus, the field at each spatial point satisfies a Langevin equation 

dV 

where the noise term is gaussian with 



dV 

3Hc/)+ — =£, (25) 



/f j 3 [ 2tt 2 C 2l Q ff5 

[Dmmt')e- Jdta ^ € i = —5{t-U). (26) 

This shows that expectation values can be obtained by averaging over different stochastic histories of the field. It 
is also evident that ([25} and (|26|) coincide with the assertions made in the introduction. Furthermore, the Langevin 
dynamics suggests that an equilibrium should be eventually reached if the potential is bounded from below as was 
first suggested in 
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To obtain the probability distribution for <f> we consider the quantity 

<t>f 



V{M<k,U) = /[^][UV]e^*[-^ +3H ^«3f)+4^"] (27) 



<t>f 



= f[Dcf,}[D y }[D^[Dp} e ^[M^y)-**(i+^v+ l Sk)-^'] 



(28) 

which represents the probability 2 to find the field value <j> at time t, given the field value <j>i at t\ in any such cell. 
Writing tp = — «p, p = —iq we have 

V fat\M = J [Dcj>][Dy\[Dp][Dp] e~ % (29) 

where the "Hamiltonian" is 

off 5 / dV\ 

H(p,q,y,4,) = ^- I p 2 -ph,Hy+ — \+q y (30) 

With p = —d y and q = — d§ and normal ordering in the Hamiltonian, the probability V will satisfy a corresponding 
"Shrodinger" equation which is nothing but the Fokker-Planck equation 



r i 9H 5 d 2 d dV d d , 

^-(^^ + 3^-,+ --^-)^ (31) 



The correlation functions generated by (|15[) can thus be computed using V '. Equilibrium is described by the stationary 
solution to (I3T1) which is easily found to be 

g 3/f 4 2 o 3/f 4 

P(tj>,y) = (32) 



with 



Therefore, any late time correlation function (£>(</>, y)) m equilibrium can then be written 



C-x - ■. ,> ! > / #dy 0(0, y) (34) 



where y = <f>. The system therefore equilibrates to a Maxwell-Boltzmann distribution. 



DISCUSSION 

The long wavelength part of a scalar field in dc Sitter, averaged over a region of physical radius R ~ -h, performs 
Brownian motion in a medium of temperature Tos = A fluctuation-dissipation relation holds between the 

amplitude of the fluctuations and the "friction" term in the Klein Gordon equation. An equilibrium is reached, 
given by (|32[) which is nothing but the Maxwell-Boltzmann distribution for the kinetic and potential energy in a 
volume ^-R%s a ^ temperature Tjjs- No thermal environment is a priori assumed and the Brownian motion emerges 



2 This is inherited from the closed time path contour in the initial path integral. 
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naturally from the effects of spacetime curvature on the quantum field. It is interesting to note that the long wavelength 
dynamics is semiclassical in that the terms suppressed by e in (fl"4| are also suppressed by powers of h 2m , as can be 
seen by reinstating h and rescaling ijj — > hip [8|. However, quantum physics affects long wavelengths through the ip 2 
fluctuation term in (1151) which induces stochastic noise. 



Although these results are in agreement with earlier stochastic treatments on the part of the distribution oc e S 1 , 
they also predict a significant average kinetic energy for the field. Equipartition is predicted for a free field, as can 

V 



also be seen from (p~9|) . while V = jcj) 4 gives 



hj?)=2{V) = ^. (35) 

The magnitude of the kinetic energy (|35|) does not of course contradict the fact that field modes are frozen when they 
enter the long wavelength sector. This kinetic energy reflects the change induced in the long wavelength field from 
the continuous influx of horizon crossing modes and is a result of choosing a fixed physical and not comoving UV 
cutoff for the long wavelength theory. It does imply however that the slow-roll assumption cannot be made for the 
long wavelength field subject to fluctuations in this case. This is a generic result of dissipative systems: no matter 
how strong the dissipation, fluctuations work to restore equipartition between potential and kinetic energy. 

A few comments are in order regarding the above results. Firstly, non-perturbative calculations in Euclidean de 
Sitter for V = \^ % E3 agree with the equal time expectation values of operators 0((j>) obtained using (|32p . 
The agreement goes beyond the leading order in powers of yf~\ if (fT3|) is expanded. This will be discussed in more 
detail in forthcoming work. Thus, it would seem that a scalar test field in Lorentzian de Sitter can be thought of 
as a non-equilibrium thermodynamic system that relaxes towards equilibrium which in turn can also be described in 
Euclidean de Sitter. A second comment relates to the infrared divergences that plague cosmological computations 



11]. We saw that in the present context such divergences are simply artifacts of perturbation theory and are removed 
if the semi-classical non-perturbative nature of the infrared sector is taken into account. Different approaches also 
suggest that such infrared divergences can indeed tamed in de Sitter and physically interpreted |12h17| |. One would 
expect that similar findings would hold for more realistic quasi-de Sitter inflationary models, also including scalar 
gravitational perturbations. These issues will be explored elsewhere. 

Before closing we would like to note that all of the above results completely ignore the backreaction of the scalar 
field on the spacetime. However, if the average energy momentum tensor obtained from (|35|) is naively plugged into 
the Friedman equations one obtains that the rate of expansion slows down according to 

H(t) = r , (36) 

leading to a decrease of the expansion rate and the temperature. Since the rate of change is small it is reasonable to 
assume that the stochastic analysis holds up even in this case. Thus, the induced stochastic kinetic energy of the field 
causes de Sitter spacetime to cool down slowly and the effective vacuum energy to decrease. For a massive particle 

this can continue only up to the point where the temperature becomes comparable to its mass after t ~ . 
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